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ABSTRACT 

The Phase-Induced Amplitude Apodization (PIAA) coronagraph is a high performance coro- 
nagraph concept able to work at small angular separation with little loss in throughput. We 
present results obtained with a laboratory PIAA system including active wavefront control. The 
system has a 94.3% throughput (excluding coating losses) and operates in air with monochromatic 
light. 

Our testbed achieved a 2.27 10 -7 raw contrast between 1.65 X/D (inner working angle of 
the coronagraph configuration tested) and 4.4 X/D (outer working angle). Through careful 
calibration, we were able to separate this residual light into a dynamic coherent component 
(turbulence, vibrations) at 4.5 1CP 8 contrast and a static incoherent component (ghosts and/or 
polarization missmatch) at 1.6 10~ 7 contrast. Pointing errors are controlled at the 1(T 3 X/D 
level using a dedicated low order wavefront sensor. 

While not sufficient for direct imaging of Earth-like planets from space, the 2.27 10 -7 raw 
contrast achieved already exceeds requirements for a ground-based Extreme Adaptive Optics 
system aimed at direct detection of more massive exoplanets. We show that over a 4hr long 
period, averaged wavefront errors have been controlled to the 3.5 10 -9 contrast level. This result 
is particularly encouraging for ground based Extreme- AO systems relying on long term stability 
and absence of static wavefront errors to recover planets much fainter than the fast boiling speckle 
halo. 

Subject headings: instrumentation: adaptive optics — techniques: high angular resolution 



1. Introduction 

An imaging system aimed at dection or charac- 
terization (spectroscopy) of exoplanets must over- 
come the large contrast beween the planet and its 
star. This is particularly challenging for Earth- 
like planets, where the contrast is « 10~ 10 in the 
visible and the angular separation is 0.1"for a sys- 
tem at lOpc. Many coronagraph concepts have 
recently been p roposed to ov e rcom e this challenge 
(see review by iGuvon et al. (|2006l )). Among the 
approaches suggested, Phase-Induced Amplitude 
Apodization (PIAA) coronagraphy is particularly 
attractive. In a PIAA coronagraph, aspheric op- 
tics (mirrors or lenses) apodize the telescope beam 
with no loss in throughput. A PIAA coronagraph 
combines high throughput, small inner working 
angle (2 X/D for 10 -10 contrast), low chromaticity 
(when mirrors are used), full 360 degree discovery 
space, and full IX/ D angular resolution. Angular 
resolution (size of the planet's PSF in the image) 
is a critical performance parameter as exoplanet 
imaging sensitivity, even if speckles have been per- 
fectly removed, is usually background-limited due 
to sky or thermal emission (near-IR and mid-IR 
imaging from the ground) or zodiacal and exozo- 
diacal light (direct imaging of Earth-like planets 
in the visible from spa ce). Th e PIAA concept, 
orginally formulated by iGuvonl §003), has since 
been stu died in depth in several sub s equent publi- 
catio n s dTraub fe Vanderbeil 120031: IGuvon et al 



20051 : JVanderbei fc Traubl 120051: iGalicher et al 
2003 : fMartinache et all 120061: Iyanderbeil""l2 006: 



Pluzhnik et al.ll2006l:lGuyon l et al.l200dlBehkov et all 
20061 : IGuvon et alJl2009t lLozi et al.ll2009h . which~ 
the reader can refer to for detailed technical infor- 
mation. 

In the first la boratory demonstrat ion of the 



PIAA concept ( Galicher et al. 2005f ). lossless 



beam apodization was demonstrated, and the field 
aberrations introduced by the PIAA optics were 
confirmed experimentally. In this first prototype, 
the PIAA acrilic optics lacked surface accuracy re- 
quired for high contrast imaging, and since this ex- 
periment did not include active wavefront control, 
the high contrast imaging potential of the tech- 
nique could not be demonstrated. In the present 
paper, we report on results obtained with a new 
system which includes reflective PIAA optics and 
wavefront control. Our prototype combines the 



main elements/ subsystems envisionned for a suc- 
cessful PIAA imaging coronagraph instrument, 
with the exception of corrective optics required to 
remove the strong off-axis aberrations introduced 
by the PIAA optics. This last subsystem has been 
designed and built for another testbed, and its 
labora tory performance is reported in a separate 
paper ( Lozi et al. 20091 ). 

The overall system architecture adopted for our 
test is presented and justified in The design of 
the main components of the coronagraphs (PIAA 
mirrors, masks) is also described in this section. 
Wavefront control and calibration are discussed in 
fJH Laboratory results are presented in fjll 

2. Laboratory system architecture 

2.1. Plate scale in a PIAA system 

Two optical conventions exist to define plate 
scale (physical distance (in meter) on the focal 
plane per unit of angle on the sky (in radian)): 
one is to follow the trajectory of the central ray in 
the optical system and measure its displacement 
on the focal plane as the source is moved on the 
sky, and the other one is to define the plate scale 
as equal to (F/Df p ) x D, where D is the telescope 
diameter and F/D f p is the beam F-ratio at the fo- 
cal plane. While the two conventions give identical 
results in conventional optical systems, they yield 
different values at the focal plane after the PIAA 
optics. This effect, due to the pupil distortion in- 
troduced by the PIAA optics, is well documented 
in previous PIAA-related publications, and leads 
to some confusion when comparing plate scale val- 
ues as both definitions have previously been used. 
In this paper, we chose to avoid using either defi- 
nition; instead, we adopt a convention where focal 
plane scale is physically defined relative to sky an- 
gle: 1 A/D in the focal plane is defined by the 
physical distance by which the PSF photocenter 
moves when the source position on the sky is off- 
set by 1 A/D. 

2.2. Coronagraph architecture 

The coron agraph architecture a dopted is a hy- 
brid PIAA (jPluzhnik et all 120061) . where beam 
apodization is shared between the aspheric PIAA 
mirrors (described in §2.4|) and a post-apodizer 
(described in §2.5p . The PIAA mirrors perform 
most of the apodization, but leave a small amount 
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x,y,z stage 




Fig. 1 . — Optical layout of the laboratory PIAA coronagraph system. The grey shaded area shows the rigid 
PIAA bench on which the two PIAA mirrors are mounted. The light source is at the upper corner of the 
figure. The focal plane mask (near the bottom, center) separates light into the imaging channel and the 
coronagraphic low order wavefront sensor (CLOWFS) channel. 



of excess light at the edge of beam (left at 0.85% of 
the surface brightness at the center of the beam), 
which is then removed by the apodizcr. Thanks 
to this hybrid approach, the PIAA mirrors are 
more easily manufacturable (less aspheric) and 
the apodizer tolerances are relaxed (the apodizer 
is not absorbing light in the bright parts of the 
beam). The hybrid design also solves t he problem 
of pr opagation-induced chromaticity (jVanderbei 
2006), which would otherwise limit contrast at 
« 10 -7 in a non-hybrid system working in a 20% 
wide band. While this second benefit was not rel- 
evant in our monochromatic experiment, it is key 
to enable high contrast direct imaging of exoplanet 
from space. 

The cost in throughput and angular resolution 
due to the apodizer are small since the apodizer 
only removes light in the fainter edges of the 
remapped beam. 

A high contrast image is formed after the 
apodizer, where starlight is blocked by the focal 



plane mask. Since the upstream PIAA optics + 
apodizer have apodized the beam with little loss in 
telescope angular resolution, the focal plane mask 
is small, with a radius ranging from approximately 
1 X/D on the sky for a 10 -6 contrast goal to ap- 
proximately 2 X/D on the sky for a 10~ 10 contrast 
goal. The focal plane mask is also part of the low 
order wavefront sensor (Guyon et al. 2009) briefly 
described in Q3.2\ which uses starlight reflected 
by the focal plane mask for accurate sensing of 
pointing errors and defocus. 

The optical layout of the laboratory experiment 
is shown in Figure [T] The light source is a sin- 
gle mode fiber fed by a HeNe laser (A = 632.58 
nm), mounted on a x,y,z stage for control of the 
input tip/tilt and focus. The PIAA system (mir- 
rors PIAA Ml and PIAA M2) creates a converg- 
ing apodized beam. PIAA M2 is chosen as the 
pupil plane for the system, and lens 11 creates a 
small image of the pupil plane onto the apodizer. 
Lens 12 reimages the pupil plane on the deformable 
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mirror, which is located in a wF/60 converging 
beam. Lenses 13 and 14 form a focal plane for the 
focal plane mask. The reimaging lens 15 is used 
to create a pupil plane and a focal plane. A Lyot 
mask is located in the pupil plane, but can be 
remotely moved out to allow the science camera 
(which is nominally in the focal plane) to move 
forward and acquire a direct pupil plane image. 
The focal plane mask reflects some of the light 
to the coronagraphic low order wavefront sensor 
(CLOWFS) came ra. A detailed descr iption of this 
device is given in iGuvon et all (|2009h . 

Once the number of lenses/mirrors and their 
relative positions was decided, a set of equations 
was coded to link together the exact location of op- 
tical elements, their focal lengths (for lenses), the 
position of pupil and focal planes, the beam size on 
each optical element, and the plate scale in the fo- 
cal planes. The equations contain all the hard con- 
straints of the experiment (for example, the DM 
must be conjugated to the pupil plane apodizer). 
Each variable was given an allowed range (for ex- 
ample beam size on the DM) or a set of allowed 
values (focal lengths of lenses constrained by what 
is available from vendors). A randm search algo- 
rithm was then used to test many possible optical 
designs and select the solutions which meet the 
criteria. This approach provided us with a flexi- 
ble tool to explore design options. The same op- 
timization code was also used to compute offsets 
in the position of several components during fine 
alignment of the system. 

2.3. Wavefront control hardware and ar- 
chitecture 

A single 32 by 32 actuators MEMS deformable 
mirror (DM) is used for wavefront control and is 
located after the PIAA optics. The wavefront con- 
trol subsystem is therefore fully decoupled from 
beam shaping effects introduced by the PIAA op- 
tics. The PIAA optics simply deliver an apodized 
beam to the wavefront control subsystem, which 
operates independently of the technique used to 
apodize the beam (remapping vs. conventional 
apodization) . While this configuration is simpler 
than a configuration where the DM is ahead of 
the PIAA optics for wavefront control, we note 
that it does not offer as wide a field of view due to 



the m agnification effect described in IGuvon et al 
<|2005h . 



Our laboratory demonstration was performed 
in monochromatic light. In this configuration, a 
single DM provides sufficient degrees of freedom to 
remove coherent diffracted light in one half of the 
field of view, regardless of how phase and ampli- 
tude aberrations in the beam are created. In a real 
coronagraphic instrument operating in broadband 
light, diffractive propagation between optics needs 
to be taken into account when designing the wave- 
front control architecture: manufacturing errors 
on optics introduce amplitude errors and wave- 
front chromaticity. In a high contrast instrument, 
such errors can only be addressed with a multiple- 
DM configuration, where the DM locations are op- 
timized to reduce residual wavefront errors over 
the spectral band used. Diffraction propagation 
effects between the aspheric PIAA optics surface 
would therefore need to be quantified when design- 
ing the wavefront control hardware. 

Two challenging aspects of wavefront control 
in PIAA coronagraphs (understanding how beam 
shaping affects wavefront control, and polychro- 
matic wavefront control) are therefore not ad- 
dressed in our experiment. 

2.4. Aspheric PIAA mirror design and 
fabrication 

The geometric remapping is performed by two 
highly aspheric mirrors, the PIAA mirrors. The 
role of the first mirror is mostly to project on 
the second mirror the desired amplitude profile, 
which is partially apodized with a faint plateau on 
the outside of the beam (see Figure [21 left). The 
apodization pr o file is described in more detail in 
Pluzhnik et al.l ( 20061 ). which includes a chromatic 
diffraction analysis of the PIAA optics used for our 
experiment. The PIAA Ml mirror acts as a con- 
verging element in the center (to concentrate more 
light in the center of the beam on PIAA M2) while 
light in the outside is diluted in a wide area of the 
beam on PIAA M2. This behavior explains the pe- 
culiar aspheric sag shown in Figure [2] on the right. 
The PIAA M2 mirror's role is to re-collimate light 
to output a beam which is apodized but free of 
phase aberrations. 

PIAA mirrors can be designed by solving a 
relatively simple differential equation when the 
input and output beams a re co llimated (the 
equation is given in iGuyonl (|2003), and also in 
Traub fc Vanderbeil ((2003) in a different form) or 
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Fig. 2. — Left: Beam apodization profile used for the design of the PIAA system. In this hybrid design, the 
apodization is not complete, and the beam surface brightness at the edge of the beam is left at 0.85% of 
the center surface brightness. This extra light will need to be removed by a conventional apodizer. Right: 
Apodization radial sag term for PIAA Ml and PIAA M2 in an on-axis configuration. The narrow region at 
the edge of PIAA Ml has a strong localized curvature, and is the most challenging feature for manufacturing 
the PIAA system optics. 



when the system is on-axis. In our laboratory ex- 
periment, the PIAA mirrors are focusing elements 
and the aspheric remapping shapes are added to 
off-axis parabolas. In this configuration, the PIAA 
mirror shapes cannot be derived from a simple dif- 
ferential equation, and they were designed by an 
iterative algorithm: 

1. Initialization: The PIAA mirror shapes are 
computed by solving the differential equa- 
tion for an on-axis system. 

2. A constant slope is added to each of the 
PIAA mirror. If there were no apodiza- 
tion, the mirrors obtained in step 1 would 
be on-axis parabolas, and this slope would 
turn them into off-axis parabolas. In the 
PIAA system, however, adding this slope 
only leads to an approximation to the off- 
axis PIAA system 

3. A 3 dimensional raytracing code is used to 
compute the beam phase and amplitude on 
the surface of a sphere centered on the out- 
put focus of the system immediately after 
reflection on PIAA M2. 

4. The difference between the measured and 
desired beam amplitude on the sphere is 



used to update PIAA Mi's shape by lin- 
ear decomposition of this residual (using 
pre-computed residuals obtained by adding 
Zernike polynomials on PIAA Mi's shape). 

5. The residual phase error measured on the 
sphere is compensated by changing PIAA 
M2's shape. 

6. Return to Step 3 with the new mirror shapes. 

This algorithm converges because changing PIAA 
M2's shape has little effect on the amplitude pro- 
file of the beam on the sphere, which is almost 
entirely a function of PIAA Mi's shape. 

The PIAA shapes were computed for a 75mm 
beam diameter at the PIAA mirrors, a 1.125m sep- 
aration from the center of PIAA Ml to the cen- 
ter of PIAA M2, and a 190mm offset between the 
PIAA Ml to PIAA M2 centerline and the input 
and output of the PIAA system (see Figure [T|) . 
In the coordinate system shown in Figure [TJ each 
mirror shape can be written as: 

z(x, y) = OAP(x, y) + sag(r) + Ejcn^ (1) 

OAP(x, y) is the off-axis parabola which would 
be the PIAA mirror shapes if no apodization was 
performed. It is a 1133 mm focal length OAP 
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with a 190 mm off-axis distance from the center 
of the optical element. This shape is identical for 
the two PIAA mirrors although the orientation is 
different. sag(r) is the apodization radial sag on 
each mirror. It is computed for an on-axis system. 
A corrective term is added to account for the fact 
that the system is off-axis (tends to for an on- 
axis system), decomposed as a sum of Zernikes 
polynomials up to radial order 7. For both PIAA 
mirrors, this correction is «95 nm RMS (excluding 
tip-tilt). 

The most challenging feature of the system is 
the small radius of curvature in the outer part of 
the beam on PIAA Ml. While the hybrid design 
adopted mitigates this problem, the radius of cur- 
vature still reaches a minimum of 155 mm near the 
edge of the mirror, at 36.1mm from the center of 
the mirror. 

The PIAA mirrors were fabricated by Axsys 
Imaging Technologies. The mirror substrates were 
initially diamond turned according to the 3-D 
prescriptions described above, and then polished 
against computer generated holograms (CGHs). 
PIAA Ml and PIAA M2 were then assembled on 
a rigid aluminum bench, aligned and permanently 
fixed to the bench. The residual system wavefront 
error was then reduced to 0.04 waves RMS by fig- 
uring PIAA M2. Two sets of PIAA mirrors (4 
mirrors total) were manufactured. 

2.5. Post-apodizer and system throughput 

In order to ease manufacturing, the PIAA op- 
tics were designed to perform most, but not all, of 
the beam apodization required for high contrast 
imaging. A more conventional apodizing scheme 
is therefore necessary to transform the beam pro- 
file at the output of the PIAA optics (solid curve 
in Figured]) into the desired beam profile (dashed 
curve in Figure [3]) . 

The post apodizer was designed in transmis- 
sion, with a series of narrow opaque rings block- 
ing light. The position and width of the rings is 
optimized to best approximate the ideal continu- 
ous apodization profile shown in Figure [3] as the 
curve labeled "Apodizer transmission" . Several 
contraints were imposed on the design to ensure 
manufacturability: no ring should be less than 0.8 
\xm wide and the gap between consecutive opaque 
rings should be no less than 5 y,m. The result- 



ing design is composed of 109 opaque rings for a 
total apodizer diameter of 3.815 mm (defined by 
the outer edge of the last opening between opaque 
rings). The apodizer was manufactured by lithog- 
raphy on a transmissive substrate. 

The post-apodizer throughput over the 3.815mm 
diameter is 96.9 %, but due to the narrow rings 
in the apodizer, some of the light transmitted is 
diffracted at large angles. The effective through- 
put of the apodizer is 94.3%, and would be equal 
to the throughput if the apodizer were continu- 
ous instead of binary. The full system throughput 
can therefore reach 94.3% (excluding losses due to 
coating) provided that the telescope pupil size on 
PIAA Ml is adjusted to the apodizer diameter. 
In practice, the telescope pupil should however 
be made slightly larger to allow for pupil center- 
ing errors, and in very high contrast applications 
(space coronagraphy) , to mitigate possible edge 
ringing effects due to Fresnel propagation. 




earn radius ( I = edge of nominal PIAA beam) 



Fig. 4. — Laboratory measurement of the apodizer 
throughput. The designed and measured radial 
throughputs are compared. 

The apodizer throughput was mesured by in- 
serting a pinhole in the PIAA output focus and 
moving the science camera in the pupil plane. The 
1 /im pinhole is used to de-apodize the beam at 
the expense of a very low throughput. Figure |4] 
shows both the measured apodizer profile and the 
designed apodization. The residual difference be- 
tween the two curves is due to the finite size of the 
pinhole (the beam before the apodizer is slightly 
apodized, so the measured profile is slightly too 
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Beam radius (1 = 37.5 nun on PIAA M2) 



Fig. 3. — Left: Apodizer design. Right: Microscope image of the outer part of the apodizer (the image 
covers approximately 1mm vertically). The outer edge of the apodizer (last transmissive ring) is at 1.9mm 
radius. The width of the individual opaque rings is 0.8/im. 



bright in the center) and the finite angular resolu- 
tion of the pupil re-imaging (the sharp edge of the 
apodizer is blurred). 

2.6. Focal plane and Lyot masks 



opaque (transmission ~ le-5) 




Fig. 5. — Microscope image of the focal plane 
mask 

The focal plane mask, shown in Figure [5j is 
used in transmission. The central part of the mask 
blocks the bright central PSF core. The radius of 
this non-transmissive central zone defines the in- 



ner working angle of the coronagraph, which is 
1.65 X/D in our experiment. A clear zone trans- 
mits the science field to the science camera. The 
shape of this clear zone of the focal plane mask 
is chosen to exclude regions of the focal plane 
where the wavefront control system cannot remove 
diffracted light. Since our experiment uses a sin- 
gle deformable mirror, diffracted light can only be 
controlled over half of the field of view. The clear 
opening in the mask is therefore D-shaped with 
an outer radius imposed by the DM actuator sam- 
pling. A slightly larger rectangular zone could also 
have been adopted, but would have imposed the 
rotation angle of the focal plane mask. The re- 
flective ring, extending from 0.8 X/D to 1.65A/Z?, 
sends some of the starlight to the low order wave- 
front sensor (LOWFS) camera. We note that the 
opaque zones of the mask are not fully opaque due 
to manufacturing considerations: their transmi- 
tion and reflection are respectively « 10 -5 and ss 
10%. 

A Lyot mask is located in the pupil plane be- 
tween the focal plane mask and the science focal 
plane. This mask is designed to block all light out- 
side the geometrical pupil and transmit all light 
within the pupil. It therefore has no effect on 
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the nominal system throughput, and its role is to 
ensure that the scattered light reaching the focal 
plane camera does not contain light outside the 
pupil. Although correcting for such light is theo- 
retically possible if it is coherent and within the 
control radius imposed by actuator sampling on 
the DM, it requires an accurate model of the coro- 
nagraph which can predict how light outside the 
pupil is affected by DM actuator positions. The 
Lyot mask was made by drilling a small hole in 
an aluminum plate, and its diameter is slightly 
smaller than the pupil size to account for align- 
ment tolerances. 

Both the focal and Lyot masks are on motorized 
stages and can be removed from the beam. 

3. Wavefront Control 

As described in section 12. 3[ wavefront control 
in our experiment is performed after the PIAA 
optics, which allows full decoupling between the 
pupil remapping introduced by the PIAA optics 
and the wavefront control algorithms. The wave- 
front calibration and wavefront control routines 
used in our experiment are therefore not specific to 
PIAA - they simply take as an input the apodized 
beam from the PIAA optics. We present in this 
section these routines as they are an essential part 
of the experiment, and the calibration techniques 
developed for this experiment are used to quantify 
the coronagraph's performance beyond the raw 
contrast achieved. Section |3~T1 describes how the 
wavefront is first flattened and the DM response is 
calibrated in a non-coronagraphic mode (no focal 
plane mask). Section [3721 briefly describes the low 
order wavefront control loop, which is detailed in 
a separate paper. In ^3.31 the main wavefront con- 
trol loop is described along with the calibrations 
used to measure the coronagraph performance be- 
yond the raw contrast. 

3.1. Initial Calibration Loop (without fo- 
cal plane mask) 

Initial calibration is performed using conven- 
tional phase diversity with no focal plane mask: 
images are acquired with the science camera in 
six positions regularly spaced from the focal plane 
to the pupil plane. An iterative Gerchberg-Saxon 
algorithm is used to reconstruct the pupil plane 
complex amplitude. As shown in Figure [6] (top 




2A/D 2 X/D 



Fig. 6. — Pupil phase (top) and focal plane image 
(bottom) when the deformable mirror is powered 
off (left) and set to its nominal position after cal- 
ibration (right). A malfunctionning actuator is 
visible on the right side of the beam. 

left), the beam quality is intially quite poor, with 
a large amount of astigmatism. The correspond- 
ing focal plane image is shown in Figure[HJ bottom 
left. 

The phase diversity routine described above is 
reapeated N times (N sa 10), with a different set 
of DM voltages applied for each phase diversity 
measurement sequence. The N phase maps ob- 
tained and the N DM voltage maps used to ob- 
tain them are then used to constrain a model of 
the DM response which consists of seven param- 
eters: Geometrical correspondance between the 
DM and the pupil image (4 parameters: x and y 
shift, scale, and rotation) and physical constants 
describing the DM behavior (3 parameters: width 
of the actuator influence function, DM dispace- 
ment for 100V applied and power index a in the 
displacement to voltage relationship with displace- 
ment cx V a ). The result of the DM calibration can 
then be used to flatten the wavefront measured 
and produce a sharp focal plane image (Figure [6j 
right). 
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Fig. 7. — Result of the automatic DM calibration. 
The DM is driven to produce a recognizable pat- 
tern on top of the flat wavefront (DM voltages, top 
left). The DM model is used to convert the volt- 
ages into a simulated beam phase (top right). The 
difference between the measured beam phase (bot- 
tom left) and the simulated beam phase is shown 
in the bottom right panel. All images except the 
DM voltage map are shown with the same linear 
phase scale. 

This initial calibration is a necessary prelimi- 
nary step for the high contrast wavefront control, 
which needs (1) a knowledge of the starting point 
(typically less than 1 radian error on the wave- 
front) and (2) a good understanding of how DM 
commands affect the pupil plane phase. The qual- 
ity of the DM calibration is shown in Figure [7] for 
a large DM offset (60 nm RMS in this example). 
The difference between the measured and simu- 
lated beam phase binned to the actuator size is 6 
nm RMS in the central 75% radius of the pupil 
(phase measurement in the outer part of the pupil 
is noisier due to lack of flux), corresponding to 
a 10% relative accuracy. The DM model relative 
accuracy is better for smaller displacements. 



3.2. Low order wavefront errors 

Low order wavefront errors are measured by re- 
flecting a portion of the bright starlight masked 
by the coronagraph focal plane into a dedicted 
camera. A detailed description of this low or- 
der wavefron t sens or (LOWFS) can be found in 
Guvon et all (|2009h . The LOWFS signal is used 
to simultaneously drive the deformable tip-tilt and 
the source position ahead of the PI A A optics. A 
key feature of the LOWFS is the ability to sepa- 
rate pointing errors (pre-PIAA tip-tilt) from post- 
Pi A A tip-tilt, which is essential to maintain high 
contrast: even a small pre-PIAA tip-tilt creates 
diffraction rings outside the IWA of the coron- 
agraph, and pre-PIAA tip-tilt errors cannot be 
compensated for by post-PIAA tip-tilt. 

When the low-order loop is closed, the mea- 
sured residual pointing error is 10~ 3 A/D, and is 
therefore small enough to be negligible in the scat- 
tered light error budget shown in ^4.31 A more de- 
tailed description of the design, calibration, con- 
trol algorithm and perfor mance of the LOWF S in 
our experiment is given in Gu von et al. ( 2009f ). 



3.3. High order wavefront control loop 

Coherent scattered light in the clear opening 
of the focal plane mask is measured by phase di- 
versity introduced on the DM. A series of focal 
plane images, each acquired with a slightly differ- 
ent DM shape, is used to reconstruct the com- 
plex amplitude and coherence of the scattered 
light. The high order wavefront control loop uses 
a linearized representation of the system in focal 
plane complex amplitude, as described in the elec- 
tric field conjugation (EFC) approach proposed by 
iGive'on et al.1 (|2007l ). The wavefront control loop 
is shown in Figure [51 and is built around the EFC 
approach. 

Prior to starting the loop, a model of the coro- 
nagraph is used to compute how each actuator 
motion affects the complex amplitude in the focal 
plane. Since this relationship is linear for small 
displacements, this model is stored as a complex 
amplitude system response matrix M resp (shown 
on the right hand side of Figure [8j of size n by m, 
where n is the number of DM actuators (ignoring 
actuators outside the pupil) and m is the number 
of pixels in the high contrast region of the focal 
plane. M resp is the linear operator which estab- 
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Fig. 8. — High order wavefront control loop, showing both the main loop and the system response matrix 
optimization loop (light shaded area). The two dark shaded boxes indicate image acquisition, which in the 
simulation mode, can be replaced with a simulated image acquisition using a model of the experiment and 
DM response. 



lishes the relationship between deformable mirror 
actuator displacements 5DM(u,v) (sampled in n 
points) and the corresponding complex amplitude 
change Sf(x, y) (sampled in m points) in the focal 
plane: 

6f(x,y)=M resp 5DM(u,v). (2) 

M resp is therefore computed by moving each ac- 
tuator of the DM in the coronagraph model and 
storing the corresponding change in focal plane 
complex amplitude in a column of M resp . 

3.3.1. Loop initialization 

For each iteration k of the loop, the first step in 
the wavefront sensing process is to acquire an im- 
age Iq with the DM shape set at the best known 
position for high contrast imaging (upper left cor- 
ner of Figured]). This first image is then used to 
choose the shapes to apply on the DM to optimize 
measurement accuracy and sensitivity. 



3.3.2. Considerations for the choice of wavefront 
sensing DM shapes 

Residual light in the coronagraphic focal plane 
is measured by adding N known wavefront errors 
on the nominal deformable mirror shape. For each 
wavefront error added, a coronagraphic image is 
taken. Intensity variations between the N images 
encode the complex amplitude of the residual light 
that needs to be measured and removed. In this 
section, we discuss how to choose the N wavefront 
errors which are sent to the DM for sensing. 

These N DM displacements are denoted 8DM^(u, v), 
with i = 1 .. N (the index k denotes the wavefront 
control loop iteration). The complex amplitude 
added to the focal plane by each of the DM dis- 
placements is given by equation [2] 

if (x, y) = M resp 6DM?(u, v). (3) 

These complex amplitude functions are referred 
to as wavefront sensing probes. The amplitude of 
the probes must be carefully chosen: if the probes 
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are too strong, the measurement is too sensitive 
to errors in the DM calibration; if they are too 
weak, the measurement is contaminated by pho- 
ton noise, readout noise and small variations in 
the incoherent scattered light. As a guideline, it is 
therefore best to choose these DM offsets so that 
the additional light (the complex amplitude focal 
plane probes) is approximately as bright as the 
light which needs to be measured. Finally, ran- 
domly modulating the probes can mitigate the ef- 
fect of calibration errors. 

3.3.3. Wavefront sensing probes 

The first probe P* is chosen to satisfy, for each 
pixel (x,y): 



\P?(x,y)\ 2 = a + a 1 I k (x,y) 



(4) 



where Jg(x,y) is the image acquired with DM 
shape DMq. If ct\ = 0, this constraint will force 
the DM shape to add a uniform coherent back- 
ground of contrast ao in the focal plane, while if 
ctQ = and a.\ = 1 it will drive the DM shape to 
add a speckle map with the same intensity as in 
the Iq image. The phase of this probe is not con- 
strained, and is chosen by an iterative scheme to 
best satisfy equation[U with SDM* as the free pa- 
rameter. We note that an exact solution to equa- 
tion|4]may not exist within the DM space, but the 
following algorithm yields a good approximate so- 
lution: 

1. Initialization: P* is computed from equation 
2] with a zero phase. 

2. Singular value decomposition (SVD) of 
M resp is used to estimate SDM* from P*. 
Eigenvalues close to zero are rejected in the 
SVD to improve stability. 

3. SDM* is clipped to avoid large DM displace- 
ments. 

4. Equation[3]is used to recompute P* from the 
clipped 6 DM*. 

5. The amplitude of P* is updated to satisfy 
equation [U but its phase is left unchanged 
from Step 4. 

6. Return to Step 2 until the iterative algo- 
rithm is stopped. 



This iterative algorithm produces simultaneously 
the probe Pf and the corresponding DM displace- 
ments. 

The second probe is chosen so that, at each 
point (x, y) in the focal plane, its amplitude is 
identical to the first probe, but its phase is off- 
set by it/2: 



P£(x,y)*iP?(x,y) 



(5) 



This 7r/2 phase offset maximizes the WFS sensi- 
tivity if the dominant source s of noises are pho - 
ton noise and readout noise (jGuvon et al. 2005). 
We note that if all DM actuators are functioning, 
there is a perfect solution to this equation, which 
can be obtained by shifting each spatial frequency 
of the DM* map by ir/2. Images acquired with 
these first two probes, together with the image Iq , 
would be sufficient to solve for wavefront errors if 
light in the focal plane is fully coherent, but at 
least one more probe is needed to unambiguously 
measure light coherence, and more probes can also 
provide the redundancy required for implementa- 
tion of the diagnostic tools described in £13.3.41 and 
g3X3 

Two additional probes have been chosen to be 



-Pf and P 4 fc 



-P 2 fc , with exact solutions 



5 DM* = -8 DM* and 5 DM* = -8 DM* respec- 
tively. 

In our laboratory experiment, we chose to also 
add 5 more probes with random uncorrelated DM 
shapes of similar amplitude than the DM displace- 
ments obtained for probes 1 to 4 above. The mean 
dispacement amplitude per actuator is measured 
on probes 1 to 4, and random displacement maps 
are produced with the same amplitude. These ad- 
ditional probes are not required for wavefront re- 
construction (probes 1 to 4 provide sufficient in- 
formation), but, as described in the next sections, 
are added to allow for calibration of M resp , mea- 
surement of incoherent light in the system, and 
measurement of coherent light variation during the 
sequence of N exposures. 

3.3.4- Solving for complex amplitude and coher- 
ence 

For each pixel (x, y) of the focal plane detec- 
tor, the complex amplitude A(x, y) of the coherent 
light leak and the intensity I(x,y) of the incoher- 
ent light leak is estimated by solving the following 
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set of equations: 

Jftz, y) = \A(x, y) + P t k (x, y)\ 2 + I(x, y) (6) 

for i — 0...N, with Pq (x, y) = (image acquired 
with nominal DM shape). 

This set of equations has three unknowns and is 
therefore overconstrained for N > 2. With N = 9 
adopted in our experiment, A(x, y) and I{x 1 y) are 
chosen to minimize : 

X 2 (x,y)= £ (Ihx,y) ~ (\A(x,y) + Pf(x,y)\ 2 

i=l..N 

(7) 

In noiseless ideal simulations, the residual x 2 
should be null. The residual x 2 m the labora- 
tory experiment includes errors due to: 

• Photon and readout noise in the I^(x,y) 
measurements 

• Variations in the light leaks during the mea- 
surement sequence. When solving for this 
set of equation, we assume A(x, y) and 
I(x, y) are static, but if they vary, \ 2 will 
increase. 

• Systematic model errors (errors in M resp ), 
which lead to errors in the estimation of the 
values of P^(x,y). If M resp is wrong, then 
the DM command sent will not produce the 
expected P*{x, y). 

For convenience, we have scaled x 2 in corona- 
graphic contrast unit. This scaling is performed 
by measuring the uncorrelated noise that would 
need to be added to A(x, y) between frames to re- 
produce the observed value of x 2 ■ In this unit, the 
observed \ 2 is approximately 4.5 10~ 8 (see 

The first contribution (photon + readout noise) 
has been computed to be a small part of the x 2 ob- 
served. We observed that increasing a% (see equa- 
tion [4]) above 1 has little effect on the residual x 2 j 
also independantly suggesting that x 2 is not due 
to detector readout noise. 

3.3.5. System response matrix optimization 

One key output of the wavefront control loop 
is the estimation of coherent light leaks (which 
should be used to compute the DM correction 
to apply for the next iteration), incoherent light 
leaks (which the DM can do nothing about) and 



the measurement of the wavefront stability during 
the measurement sequence. All these quantities 
depend upon a reliable estimation of M resp . It 
is for example possible, if M resp is wrong, to ob- 
tain a low value of the residual coherent light and 
think the system has converged to a good contrast 
value, while in fact a significant amount of coher- 
ent light remains. This last issue is mitigated, 
but not entirely addressed, by continuously vary- 
ing the probes (as this error is a function of the 

I(x iR^S^° eS cnosen ' anc ^ wn ^ average to zero in the lin- 
' ear approximation used in this work if the probes 
are randomly chosen). An error in M resp would 
first appear as a large x 2 value for the solution of 
equations [5] 

To address this, we have added a M resp opti- 
mization loop within our control loop. For each 
iteration k, the derivative of: 



X 



(8) 



■ivy 



with each element of M resp is computed (this is 
a total of 2 x n x m derivatives, as the derivative 
is computed for the real and imaginary parts of 
each element of the M resp matrix) . This derivative 
is computed from equation [7] by replacing by 
M resp 8DM]° (equation El). 

At each iteration of the wavefront control loop, 
M resp is then slightly modified in order to reduce 
X 2 : for the 10% of M resp elements showing the 
largest derivative against x 2 > the value of the ma- 
trix element is moved in the direction indicated 
by the derivative by 0.1% of the RMS value of ele- 
ments in the corresponding column of M resp . This 
algorithm was first tested on simulated data with 
an initial M resp estimate which was different from 
the actual M resp used in the simulation for com- 
puting the images. This test showed that M resp 
did converge toward the true M resp , and that the 
X 2 value decreases as a result. Convergence is very 
slow, due to the large number of coefficients in the 
M resp matrix, requiring several hundred iterations 
before a significant improvement in x 2 is observed. 



3.3.6. Correctic 



ed to the DM 



The li near electric field co njugation (EFC) al- 
gorithm (jGive'on et al.ll2007l ) is used to cancel co- 
herent scattered light. This algorithm uses the lin- 
earized coronagraph model which is also used for 



12 




the measurement step described above. The sys- 
tem response matrix is inverted to build a control 
matrix which is multiplied to the coherent light es- 
timate to yield the DM shape offset to be applied. 



Regul arization schemes proposed bv lGive'on et al 



(|2007l ) were used to improve the loop stability and 
convergence speed. 

4. Laboratory results 

4.1. Alignment and Apodization measure- 
ment in the pupil plane 

The pupil plane apodization map is measured 
by placing the science camera in the pupil plane. 
In our experiment, the pupil is conjugated simul- 
taneously to the PIAA M2 mirror, the apodizer, 
and the DM. Alignment is necessary to ensure that 
these three planes are conjugated and that their 
relative scales are correct. The camera positions 
for which conjugation to these planes is achieved 
are measured and the corresponding pupil scales 
are derived from the images. These six numbers 
are then fed to an optimization routine which com- 
putes the offsets to be applied to all movable op- 
tical elements after PIAA M2 to meet the conju- 
gation and scale requirements. A few iterations of 
this sequence were sufficient to converge. 

In the fine alignment step, the apodizer alone 
is moved. The pupil image is compared to a sim- 
ulated pupil image where the relative scale and 
lateral offset between the PIAA apodization and 
the apodizer transmission map are free parame- 
ters. The values of these three parameters which 
give the smallest residual difference is then used to 
guide fine alignment of the apodizer. Fine tuning 
of the scale between the apodizer and the PIAA 
apodization is possible because the beam at the 
apodizer is non collimated: apodizer motion along 
the optical axis changes this scale. Figure [9] shows 
both the pupil apodization profile measured after 
alignment and the theoretical profile. 

4.2. Imaging with a non-corrected PIAA 
system 

Figure |TD] shows the system on-axis PSF in 
imaging mode (no coronagraph focal plane mask). 
The on-axis PSF is similar to an Airy function 
without the Airy rings beyond 1.22 X/D. 

While the on-axis image is sharp and exibits 




Beam radius ( I - edge of die beam) 



Fig. 9. — Measured apodization profile, compared 
the to the beam profile the experiment was de- 
signed to deliver. A 2-D image of the beam is also 
shown. 



high contrast, our laboratory system did not in- 
clude the inverse PIAA optics necessary to cor- 
rect for the strong off-axis aberrations introduced 
by the forward PIAA optics. These inverse optics 
do not need to be of coronagraphic quality, and 
can be a small set of lenses. A laboratory demon- 
stration of wide field correction with inverse PIAA 



optic s is described in a separate paper (jLozi et al 
20091) . 



Without inverse PIAA optics, the image of an 
off- axis source rapidly changes shape as the source 
moves away from the optical axis. Figure [11] shows 
the measured off-axis PSFs and the result of a 
simulations u sing a remapping of the beam phase 
( Guvonl 2003^ 1 . With such a strong field aberra- 
tion, measuring the focal plane plate scale is chal- 
lenging and its value is a function of the metric 
used. In this paper, we choose to adopt the non 
coronagraphic PSF photocenter to measure plate 
scale. In FigurefTUl the 4 X/D ring therefore shows 
where the photocenter of the PSF would be if the 
source was 4 X/D from the optical axis of the en- 
trance telescope. The bright PSF core at this sep- 
aration would be slightly outside the ring, but the 
fainter assymetric diffraction arcs of the off-axis 
PSF would be inside the ring. The same photo- 
center metric is used to measure the angular sizes 
on the focal plane mask given in H2.6I The plate 
scale on the coronagraph camera was measured 
by 2-D fitting of measured and simulated off-axis 
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Fig. 10. — On-axis PSF without focal plane mask 
(log scale indicated at the bottom of the figure). 
The vertical bleeding feature extending downward 
of the PSF core is a camera artifact, and is re- 
moved when the focal plane mask is used. A faint 
ghost due to the entrance window of the camera 
is visible just beyond 2 X/D. 

PSFs at several angular separations. 

4.3. Coronagraphic performance: mea- 
surements 

Coronagraphic performance is measured with 
the focal plane and Lyot masks in the beam and 
both the LOWFS loop and high order wavefront 
control loop closed. Figure fTS] (upper left) shows 
a raw image from the science camera. Most of the 
light reaching the camera is due to partial trans- 
mittance (at the 10~ 5 level) of the core of the focal 
plane mask, which produces a central peak in the 
image. The clear 1.65A/D to 4AX/D opening in 
the focal plane mask is visible below this central 
peak. 

Residual light is decomposed in two compo- 
nents by the wavefront sensing sytem, using the 
approach described in £)3.3I 

• An incoherent component composed of light 
which does not interfere with light extracted 
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Fig. 11. — Off-axis PSF without focal plane mask, 
for a 1.7 X/D off-axis distance. A cross indicates 
the optical axis in the simulated image (left). Both 
images are shown to the same brightness scale 
(bottom). 

directly from the central PSF core. This 
component appears to be mostly stable in 
structure and varies from iteration to iter- 
ation between 1.5 1CP 7 and 2 1CP 7 in con- 
trast. 

• A coherent component which is used to com- 
pute the DM shape for the next iteration. 
This component is at approximately 510 -8 
contrast and is decorrelated on timescales 
above the response time of the wavefront 
control loop. The estimation of this com- 
ponent varies greatly from iteration to itera- 
tion, and is sometimes below 10~ 8 contrast. 

We also measured the \ 2 of the solution, and found 
it to correspond to a change in coherent light at the 
4.5 10 -8 contrast level during the measurement se- 
quence, which is 10s long. The M resp optimization 
loop described in ^3.3.5l eave no noticable improve- 
ment in x 2 , suggesting that the inital calibration 
led to a good estimate of M resp and that the ob- 
served x 2 is indeed dominated by fluctuations in 
coherent light at the 4.5 10~ 8 contrast level. 

Visual inspection of the coherent and incoher- 
ent portions of the light strongly suggests that this 
decomposition was successful. First, the images 
obtained show very little high spatial frequency 
noise. Although the reconstruction is performed 
independantly for each pixel, both the coherent 
and incoherent components show speckles and fea- 
tures covering several pixels in these oversampled 
images. Second, the ghost to the lower left of the 
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Coherent portion 
Average Contrast = 4.48e-8 




2 X/D 



Incoherent portion 
Average Contrast = 1.63e-7 




2 X/D 



Complex Amplitude average 
of coherent portion over 1300 
loop iterations 
Average Contrast = 3.5e-9 
(Contrast scale xlO in image) 




2 X/D 




Fig. 12. — Top left: Raw coronagraphic image. A decomposition of the scattered light into a coherent 
component (top right) and incoherent component (bottom left) shows that the raw contrast is dominated 
by incoherent light. The coherent bias over a long period of time is shown in the bottom right. 
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optical axis was properly analysed by our algo- 
rithm. This ghost, also visible in figure \W\ is due 
to internal reflection in the window of the CCD 
camera. Because the two surfaces of the window 
are not perfectly parallel, interference fringes can 
be seen across this ghost. Although the fringes 
are coherent, they are too thighly spaced for the 
wavefront control system to remove. They are the 
only known feature in the image that (1) should 
be seen in the coherent image, (2) should not be 
seen in the incoherent image, and (3) cannot be re- 
moved by the wavefront control system. As shown 
in figure IT2l our coherent /incoherent analysis cor- 
rectly identifies these fringe as coherent, and only 
a minute fraction of the fringes is present in the in- 
coherent portion of the image (likely due to small 
motion of the fringes during the measurement se- 
quence) . 

4.4. Analysis 

The analysis results are summarized in table [1] 
and discussed in this section. 

4--4-1. Incoherent light component 

The raw contrast is dominated by a very sta- 
ble incoherent component at 1.63 1CP 7 , which 
is likely due to ghosts and/or polarization miss- 
match. Given the high number of air-glass sur- 
faces (twenty) , including some which are not anti- 
reflection coated, a ghost at this level would not 
be too surprising. A possible source of polarization 
missmatch is the ring apodizer which has many 0.8 
/im wide annular rings. The variations observed in 
the estimate of the incoherent light are due to the 
coherent light variation during the measurement, 
which affects the incoherent estimate. The data 
obtained is compatible with a fully static incoher- 
ent background, as would be expected from ghosts 
and polarization effects. We note that a fast vary- 
ing coherent component could also be responsible 
for this light if it is varying sufficiently rapidly to 
appear incoherent within « 1 second. 

4.^.2. Coherent light component 

During the measurement sequence, coherent 
light varies by w 4.5 10 -8 in contrast due to turbu- 
lence or vibrations in the system, as shown by the 
X 2 analysis. The coherent light leak estimate over 
a 10 second period is « 5 10~ 8 , which is at the level 



expected from the 4.5 10 -8 variations shown by 
the x 2 analysis. The large variation, from iteration 
to iteration, observed in the coherent light residual 
is due to the turbulence/ vibrations in the system. 
We note that the lucky iterations where the 10 
second averaged coherent light is estimated below 
10~ 8 are artefacts of the time averaging during 
the measurement period: even during these lucky 
periods the coherent light did vary by w 4.510 -8 
between the 10 frames of the measurement. 

4-4-3. Static wavefront aberrations 

The wavefront control loop successfully removes 
static coherent speckles. Over a 4hr period of 
time, we have measured the static coherent speck- 
les to be at the 3.5 10 -9 contrast level. Except for 
a known ghost on the camera window, we could 
not detect any residual bias above this level in the 
residual time- averaged coherent light. 

The quality of the DM calibration (better than 
10% for small displacements) and the use of 10 
probes chosen to avoid systematic bias ensures 
that within a single wavefront measurement, the 
coherent wavefront error is below 0.1 x \/To ~ 
0.03 times the DM probe amplitude. Since DM 
probes were chosen to be at the 10~ 7 contrast 
level, the statistical measurement error due to 
DM calibration errors is expected to be below 
the 10~ 7 x 0.03 2 « 10~ 10 contrast level for each 
wavefront measurement. Decorrclation between 
wavefront measurements further reduces the effect 
of DM calibration errors on the static wavefront 
aberration measurement accuracy. DM calibra- 
tion and other potential sources of errors are fur- 
ther reduced by the system response matrix opti- 
mization described in £13.3.51 At the 3.5 10 -9 con- 
trast level, the measurement of static wavefront 
aberration is therefore robust. 

5. Conclusion 

The results obtained in this experiment are es- 
pecially encouraging for ground-based coronagra- 
phy. The 2 10~ 7 raw contrast we have achieved 
in the 1.65 to 4.4 X/D angular separation range 
already exceeds by two orders of magnitudes the 
raw contrast that can be hoped for i n even a the- 



oretic ally ideal Extreme- AO system (jGuvon et al 
120051 ). We note that with a more careful opti- 
cal design and anti-reflection coated optics, our 
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Table 1 

Contrast budget (Averaged over 1.6 to 4.4 X/D) 



Term 



Likely origin 



Value 



Calibration 



Incoherent light, static 
Coherent, variable in t < 10s 
Coherent residual (in 10s) 

Coherent, static (in 4hr) 



optical ghost, polarization 
turbulence, vibrations 



uncorrectable coherent light 



1.63 10 w 
4 10~ 8 
5 10" 8 (typical) 
7 10" 9 (best) 
<3.5 10" 9 



Incoherent portion of WFS data 
residual from WFS reconstruction 
Coherent portion of WFS data 

time-averaged coherent light 



experiment could probably have reached 5 10 
raw contrast (level of coherent light leak in the 
current experiment). More importantly, we have 
demonstrated that with the coronagraph + wave- 
front control architecture adopted in our experi- 
ment, static coherent speckles can be pushed down 
very low (3.5 10 -9 ) in long exposures. Our system 
successfully removed long term correlations in the 
coherent speckles, and their averaged level in long 
exposure was reduced with a 1/yT law. The com- 
bination of a high performance PIAA coronagraph 
and a focal-plane based wavefront control there- 
fore appears extremely attractive for ground-based 
Extreme-AO. In that regard, our experiment has 
been a successful validation of the key technolo- 
gies and control algorithms of the Subaru Corona- 
graphic Extreme- AO (SCExAO) system currently 
in assembly. The major differences between the 
SCExAO PIAA coronagraph and our laboratory 
prototype are (1) the need to design and oper- 
ate a PIAA coronagraph on a centrally obscured 
pupil with thick spider vanes and (2) the need for 
corrective optics to recover a wide field of view. 
These two requirements have been validated in 
a separate laboratory experi ment using the fi nal 
SCExAO coronagraph optics ( Lozi et al. 20091 ). 



Our experiment was limited at the 2 10~ 7 con- 
trast by an optical ghost and at the 5 10~ 8 con- 
trast by turbulence or vibrations. The PIAA coro- 
nagraph could therefore not be tested to the con- 
trast level required for direct imaging of Earth- 
like planets from space (approximately 10 -10 ), al- 
though several key concepts were demonstrated, 
including simultaneous operation of a low-order 
wavefront sensor using starlight in the PSF core 
and high-order wavefront sensor using scattered 
light in the science focal plane. New calibration 
schemes which will be very useful for high contrast 



coronagraphy were also developed and validated, 
such as the system response matrix optimization 
loop, which can executed as a background task to 
fine-tune the system. 

PIAA coronagraph technologies for high con- 
trast space applications are now being actively de- 
veloped and tested at NASA Ames Research Cen- 
ter and NASA Jet Propulsion Laboratory. A new 
set of PIAA mirrors was recently manufactured 
to higher surface accuracy than the ones used for 
this demonstration, and is being integrated within 
the High Contrast Imaging Testbed (HCIT) vac- 
uum chamber at NASA Jet Propulsion Labora- 
tory. We note that the HCIT chamber has already 
demonstrated stability to the 10~ 10 contrast wit h 
a Lyot-type coronagraph ( Trauger k. Traubll2007 ). 
The experiment described in this paper served as 
a precursor to this new step, which is aimed at 
reaching higher contrast (minimum goal of 10~ 9 ) 
in broadband light using a two-DM wavefront cor- 
rection. In parallel to this effort, a highly flexible 
high stability air testbed at NASA Ames Research 
Center is coming online to further explore technol- 
ogy and architecture trades for PIAA systems. 

In addition to pushing the contrast further, fu- 
ture laboratory demonstration of the PIAA coro- 
nagraph will explore chromaticity and wavefront 
control issues unique to a PIAA coronagraph ar- 
chitecture. The monochromatic experiment de- 
scribed in this paper did not address these impor- 
tant points, and should therefore be considered as 
only a first step toward validation of the PIAA 
coronagraph technique for high contrast space- 
based exoplanet imaging mission. 

This research was conducted with funding from 
NASA JPL and the National Astronomical Ob- 
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servatory of Japan. Technical input and advice 
from the members of NASA Jet Propulsion Lab- 
oratory's High Contrast Imaging Testbed (HCIT) 
team, NASA Ames Research Center's coronaraph 
team, and Princeton University's coronagraph 
team have been of considerable help to conduct 
this work, both for design/simulations and labo- 
ratory implementation. In addition to providing 
laboratory space and infrastructure, Subaru Tele- 
scope made this research possible through major 
contributions from its technical staff (electronics, 
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